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Abstract 

We study the phase transition in the honeycomb dimer model (equivalently, mono- 
tone non-intersecting lattice path model) . At the critical point the system has a strong 
long-range dependence; in particular, periodic boundary conditions give rise to a "res- 
onance" phenomenon, where the partition function and other properties of the system 
depend sensitively on the shape of the domain. 

1 Introduction 

We study a model of monotone non- intersecting lattice paths in 1? . Applications of this 
model include random surfaces [2, 23, 15], magnetic flux lines in superconductors [32, 3], 
and a number of other physical phenomena [20, 11, 8], including spin-domain boundaries of 
the three-dimensional Ising model at zero temperature [4]. (See also [22, 13].) 

Let -Rm,,n be a domain consisting of the m x n rectangle in with periodic boundary 
conditions {Rm,n is a graph on a torus). On Rm,n consider configurations consisting of 
collections of vertex-disjoint, monotone northeast-going lattice paths (or rather loops, since 
the paths are required to eventually close up, possibly after winding several times around 
the torus: there are no "free ends"). See Figure 1. We do not restrict the number of disjoint 
loops but rather give a configuration an energy EbNb + EcNc where Nh is the total number 
of "east" steps of the paths and Nc is the total number of "north" steps of the paths. The 
Boltzmann measure fi at temperature T is the probability measure assigning a configuration 
a probability proportional to e-i^b'^b+EcNc)/T ^ study these Boltzmann measures near the 
critical temperature T, which is the temperature at which e~^*'^'^ + e~^''^'^ = 1 [14]. Letting 
b = e~^^/'^ and c = e~^''^'^ ^ a configuration has a probability proportional to h^^c^". 

This process is equivalent to another well-known model, the model of dimers (weighted 
perfect matchings) on the honeycomb lattice H [30, 31]. Weight the edges of the honeycomb 
lattice a = 1, 6, or c according to their direction as in Figure 1. See Figure 1 for an illustration 
of the weight-preserving bijection between dimer configurations and lattice paths. 

The research leading to this article was conducted in part while the first author was visiting Microsoft. 
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Figure 1: The classical bijection [30] between perfect matchings of the hexagonal lattice and 
non-intersecting north-east lattice paths, shown foramxn = 7x6 region. Edges of type 

'a' in the hexagonal lattice Hm,n are contracted to obtain -Rm,n- The dimer configuration has 
weight a^^b^'^c^'^, as does the lattice path configuration (when a = 1 or else one introduces a 
factor of a for each vertex not in a lattice path) . 



Let Hjn,n be a finite graph which is a quotient of H by the horizontal and vertical 
translations of length m, n respectively as in Figure 1. Let fj,m,n — IJ'm,nici',b,c) be the 
Boltzmann measure on dimer configurations on the toroidal graph Hm,n- This measure 
assigns a dimer configuration a probability proportional to a^'^b^^'c^" , where there are 
edges of weight a, Nj, edges of weight 6, and Nc edges of weight c. The exact probability 
is a^'^h^^c^" / Z, where the normahzing constant Z is called the partition function of the 
system. 

As m, n — > oo the measures fim,n{(i, b, c) have a unique limiting Gibbs measure /x(a, b, c) 
[25, 6]. The measure fi{a,b,c) is a measure on dimer configurations on H, and is well- 
understood for all a, b, c [6]. As a, b, c vary the measure /x(a, b, c) undergoes a phase transition 
-the "solid-liquid" transition- when one of a, b, or c is equal to the sum of the other two, 
for example a — b + c. When a> b + c the system is frozen: with probability 1 only edges of 
weight a are present. See Figure 2. When a, b, c satisfy the strict triangle inequality (each is 
less than the sum of the others) dimers of all types are present in a typical configuration of 
/x(a, b, c). Likewise in the lattice path model, when a >b+c there are no lattice paths. When 
a, b, c satisfy the strict triangle inequality, the lattice paths in a configuration are dense, that 
is, on average they lie within a constant distance of one another. 

The finite- volume measures yUm,n(a, b, c) are less well understood near a — b + c. When 
a = b + c the lattice paths exist but are spread out; as we shall see, the average distance 
between strands is on the order of the square root of the system size. The ratio n/m 
imposes fairly rigid entropic constraints on the way these loops can join up. Surprisingly 
these constraints become stronger with increasing system size, so that the scaling limit has 
nontrivial structure. 

Our main result is a computation of the partition function as a function of a, b, c, m, n, 
for parameters near a = b + c. The partition function Z (normalized by (area)^/^) is largest 
when nb/ (mc) is rational with small numerator and denominator. 

For example we have 
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Figure 2: The partition function is shown in the upper panels, and the expected number of 
edges of type 'c' is shown in the lower panels. The solid and liquid phases are shown in the 
left panels, and the melting transition between between these phases is shown in the right 
panels. The formulas on the left are derived in [14, 29, 6] and are not used here. The formulas 
on the right follow from Theorem 2. The melting transition depends quite sensitively on the 
aspect ratio of the region (see Figure 3 and Theorem 1). 
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Theorem 1 When a — l,b — c — 1/2 and n/m — p/q in lowest terms we have 

as m and n tend to oo while p and q remain fixed (here ( is the Riemann zeta function). 
Furthermore the number Nf, of 'c '-type edges tends to a Gaussian with expectation 

(mn)»/\(l/2)(l-2V^) 

(Both C(l/2) and C(-l/2) are negative.) 

We see that the system greatly prefers domains of rational modulus over those with 
irrational modulus. Here by rational we mean, a lattice path with average slope c/b will 
close up after winding a small number of times around the torus. In such a case the number 
of loops can be large, whereas in an irrational case each loop must wind many times around 
before closing up (unless it pays a large entropic cost). 

One can think of the partition function Z, taken as a function of m and n, as an indicator 
of rationality of n/m. See Figure 3 which plots logZ/(area)^/^ and (A^c)/(area)^/^ as a 
function of log(n/m) near area = 10^ for a = 1 and 6 = c = |. 

This is not the only interesting phenomenon in this model. For a rational domain, there is 
a non-trivial behavior as we vary a, b, c away from the critical point. Letting A = (c/ (a — 6))", 
the partition function as a function of A has an infinite number of nonanalyticities (in the 
large-n limit) which correspond to abrupt changes in the winding number of curves in a 
typical configuration. That is, the curves "ratchet" at well-defined values as A increases: see 
Figure 5 and § 5.4. We did not prove (although we believe) that at a typical value of A the 
curves are in a well-defined integer homology class (i.e., have well-defined winding number 
around the torus), and this homology class changes at discrete values of A. We prove only 
that the Z2 homology class (winding number modulo 2) changes at these well-defined points. 
We also show that, away from these transition points, the number of edges is a Gaussian; at 
the transition it is a mixture of two Gaussians, coming from a Gaussian for each homology 
class. 

It would be very interesting to study this same model on higher-genus surfaces. On 
higher-genus surfaces (with translation structures having conical singularities) it would be 
very useful to be able to detect "rationality" , in the form of bands of parallel closed geodesies: 
this is an important problem in billiards [26] . Moreover the ratcheting phenomenon must be 
significantly richer in the presence of a non-abelian fundamental group. 

2 Review of log Z and the distribution of Nc 

The partition function Z is 



Z^ J2 a^-'fe^^c^^ 

configurations 
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Figure 3: Resonant spikes in logZ (upper panel) and in the edge density (lower panel) for 
very large regions when a — \ and 6 = c = 1/2. The aspect ratio is plotted on a log scale. 
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It contains all the information about the distribution of the total number of edges of each 
type. Wc can view Z as a polynomial in c, where the coefficient of c^" is the weighted sum 
of configurations which contain Nc edges of type 'c'. Thus the expected number of edges of 
type 'c' is {Nc) = f Similarly, the £th moment of the number of edges of type 'c' is given 

by 



The ith cumulant Kg, is defined by 



d d 

^i^^a ^' 

oc oc 

which we will in effect estimate later. Here we recall several properties of cumulants [7]. Of 
course {Nl) — Ki. The variance in the number of edges of type 'c' is K2: 

d dlogZ _ d cdZ 

dc dc dcZ do 

_ c d dZ c dZ dZ 

Z dc dc Z'^ dc dc 

= (iVc) - (Nc)' = a\Nc). 

Later we will use higher moments to show that the number of 'c'-type edges tends to a 
Gaussian, and to this end we express these moments in terms of the i^^'s. Since 

when (iVf) is expressed as a polynomial in the variables f^i, . . . , Kg, we may calculate {N^'^^) 
from (N^) by replacing each monomial Ki-^Ki^ ■ ■ ■ Ki^ of {N^) with 

KiKi^Ki^ ■ ■ ■ Ki^, + Ki^+iKi^ ■ ■ ■ Ki^ + Ki^K^^^i ■ ■ ■ Ki^, + ■ — h K^^Ki^ ■ ■ ■ Ki^+i. 

Thus for example we have 

(Nl) =K, 
{N'^) =K, + Kl 

{Nl) =K^ + ?>K2Ki + Kl (1) 
{Nl) =K4 + AKsKi + ml + &K2KI + Kt 

{Nl) =K^ + hK^K^ + l<dK^K2 + im^Kl + IhKlK^ + l^dK^Kl + Kl 

These polynomials {Ni) — Yj{Ki, K2, ■ ■ ■) are the complete Bell polynomials [7]. We see 
that (A*"^) contains a monomial for each partition of £, and the coefficient associated with 
partition with distinct part sizes Si > S2 > • • • > Sfc and parts of size Sj is 
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It will be more useful to work with moments about the mean rather than moments about 
the origin. Note that if we replace Z with Z* = Zc~'^, then the above derivation shows us 
how to express {{Nc — /lY) in terms of the K^^s defined by 



As — Ki — fj, and — Ke ior i > 1, upon substituting // = (Nc) — Ki we see that the 
ith moment of Nc about the mean may be obtained from the above expressions for (N^) by 
deleting all monomials that contain the variable Ki. 



3 Product form of the partition function 

We compute an expression for the partition function as a function of a,b, and c. We are 
interested in approximating Z to within l+o(l) multiplicative errors when a — 1, b,c & (0, 1) 
and 6 + c is close to a. The interesting range is when (c/ (a — 6))" is of constant order, that 
is, a — b — c = 0{l/n). We define A — {c/{a — b))'^ as the natural parameter measuring 
proximity to the critical point. 

In what follows we always set a = 1, although we keep using a for notational convenience. 

Recall Hjn,n: the m x n hexagonal toroidal graph shown in Figure 1. By Kasteleyn [14] 
(see also [27, 12, 24, 18, 19] for extensions and further developments), the partition function 
Z — Z{a, b, c) for dimer coverings of is a sum of four expressions, 

Z=^(-Zoo + Zoi + Zio + Zn), (3) 

where Z^t is the determinant of a signed version of the adjacency matrix of Hm,n, and counts 
dimer coverings with a sign according to the homology class (in i^i (torus, Z2) = Z|) of the 

corresponding system of loops, as follows. Let N{ex,By) denote the total weight of dimer 
coverings whose corresponding loops have mod 2 crossings of the line x = and Sy mod 2 
crossings of the line y = 0. Each Z^^ is a linear combination of the N{ex, £y) with coefficients 
±1 as follows: 

iV(0,0) iV(l,0) 7V(0,1) 7V(1,1) 



Zoo 


+1 


-1 


-1 


-1 


Zio 


+1 


+1 


-1 


+1 


Zqi 


+1 


-1 


+1 


+1 


Zu 


+1 


+1 


+1 


-1 



(4) 



Note three important facts, which follow from this table: 

Proposition 1 The sum of any two of —Zoo-, Zoi, Zio, Zu has only nonnegative coefficients. 
The difference between any two of —Zoo, Zoi, Zio, Zu is bounded by the sum of the other two. 
The difference between the coefficients ofa°'b^c'^ in any two of —Zqo, Zqi, Zio, Zu is bounded 
by the sum of the coefficients of a°'b^c'^ in the other two. 
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Kasteleyn [14] evaluated the determinants Z^^ by multiplying eigenvalues obtained through 
Fourier analysis, giving 

Z^T = JJ^ [a + 62; + cw] 

(-2)'" = (-l)'' 

{-wr={-\y 



[(a + hzY - {-\yd' 

n (a + n 



1-f-ir 



a + bz 

c 



a + bz 



(5) 



n 

We can ignore the 6™ term which is exponentially smaller than a™. When & + c is close to 
a, unless z is close to —1, |a + 62;| will be greater than c; in particular |a + 62;|" is exponentially 
larger than c". So we can ignore the factors in the products for which z is not close to — 1. 
We will expand the remaining factors near z = —1. As 2; is a root of unity, let z — Zk — — 6*^*=, 
where 9k = 2Txk/m for k G + (j/2. Of course /c = /c + m, so when doing series expansions 
we can take — m/2 < k < m/2. 

Define r^, 4>k by 



1± 



1 ± 



a — 6 



^ a + bzk , 
so that 

T},^ {a-bY\a + bz\^y and 0jk = arg(a + 62;^)"". 
We make the simplifying assumptions 6/a = ©(l), 1 — 6/a = ©(l), and n = ©(m). Then 

+ O(n^fc) ) = exp(-eA;2 + 0(A;Vm3)) 



Tfc = exp I —nQ^ 



''2{a-bY 



where we have defined 
and 

where we define 

Letting A — (^)", we have 



lii^nab 
ra^ia — by 



0(l/m), 



nO^b 
a — b 



+ nO{9l) = (l)k + 0{kym^) 



= 



27rnb 
m{a — b) 



In logarithmic form we can write 

logZ,, = - V Lii {{-lYAvke'^") +n\og{a^ - {-lyb^) 



(6) 



(7) 



negligible 
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where the polylogarithm function Lij, is defined by Li,^(z) = X^^i ^ for jz] < 1 and by 
analytic continuation elsewhere (see Appendix A for background on polylogarithms) . Here 
the second term is essentially zero, and the terms in the summation become negligible when 
I A; I is larger than 0(v^). 

In this expression for log Zfj^ we need only keep track of its real part: from (6) we see 
that (since b < a) Z^^t is real and nonnegative if A < 1, and real and nonnegative if ^4 > 1 
except for Zqq which is strictly negative. In particular when A > 1 the expression (3) is a 
sum of nonnegative terms. When A < 1 we shall see that Zqo is negligible compared to Z 
so its sign is irrelevant. 

More generally we find, using (7), z-^ l^iuiz) = Li^-i{z), and {cd/dc)A — nA, that 



(8) 

(+nlog(a"^-(-l)'^6™) if£ = 0). 

" V ' 

negligible 



'dc-'-'dc^""^^^-^^ 



4 Rational tori 

The expressions (8) are non-trivial to evaluate, mostly because they are describing behavior 
which depends sensitively on the parameters defining the system. In this section we compute 
the asymptotics of (8) for "nearly rational" domains. 

We consider the toroidal hexagonal graph i/m.n to be "nearly rational" when 0/(27r) = 
^^^^^^ is close to a simple rational p/g, where p and q arc relatively prime integers. (Note 
that this depends not only on m, n but also on a, 6, c.) We keep p and q fixed as m and n 
tend to infinity, and by "close" we mean that nhj (mc) —p/q is not large compared to 1/y^qn. 
We introduce the parameter a to measure the closeness of nh/{mc) to p/q: we define a by 

where W = ^/q^/{7^p) = 0(^^)- The interesting range is when a is of constant order. 



We will determine the asymptotic shape of the resonant peaks in Figure 3 as functions 
of a and A. 



4.1 Spokes, spirals, and clouds 

For fixed p and q, as the area mn gets large, the terms 1 ± ArfeC*'^*' 1 ± ArkC^^^^^'^^'' 
accumulate on q different spokes (or radii) of a circle with radius A centered at 1. If 0/(27r) 
is only approximately p/q, then each of these spokes becomes a spiral, which spirals out from 
1 when k is negative and increasing, and then spirals back in towards 1 when k is positive 
and increasing. If 0/(27r) is far from a simple rational, then the terms 1 ± Ar^e*'^'' form a 
cloud within the disk of radius A centered at 1. In our analysis for nearly rational tori, we 
will assume that 0/(27r) is sufficiently close to a simple enough rational p/q that the terms 
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Zoo for 60000 x 40000 



Zoo for 60000 x 40500 



Zoo for 75025 x 46368 



Zii for 60000 X 40000 



Zii for 60000 X 40500 



Zii for 75025 x 46368 



Figure 4: The multiplicands for Zoo and Zn for a domain whose aspect ratio is (a) sim- 
ple rational (b) nearly simple rational (c) far from simple rational. The multiplicands are 
complex and accumulate towards 1. 



1 ± Arkc"^'' form what appear to be q continuous spokes or spirals (in a sense we define more 
precisely below). See Figure 4. 

It is useful to re-express (8) to reflect the presence of the q spirals. Since ^ 27r|, every 
qth term lies on a given spiral, so we break the sum apart into q different sums, one for each 
spiral. Of course q may not evenly divide m. The most convenient way to re-express the 
summation is 



E 



keZm+cr/2 



Since spirals arc continuous objects rather than discrete sets of points, we wish to approxi- 
mate XlueZm with Jp™ du. To this end we subtract 2'Kpu from (pj+qu, then the angle ( mod 2-n) 
is unchanged for integer m, while for continuous u it is slowly varying so we can hope that the 
integral approximates the sum. Thus we re-express (8) as (ignoring the negligible 0(n6™) 
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error term when £ = 0) 



jeZq+o-/2 ueZm 

Later we will quantify the error introduced by approximating these sums with integrals, and 
show it to be insignificant so long as the points appear to line up on q spirals which miss the 
singularity. For now we proceed with the integral approximation and simplify it: 

changing variables k — j ^ qu and using the fact that the integrand is periodic, 

1 



r Jo 



qm 



using the replication formula ^ ^^^9=1 ljij,{ujz) = hiy^z'^) and the fact that gcd(p, = 1, 

? Jo 



and again using the periodicity of 4>k and r^, 

(nnY 

^ / Lii_^ ((_i)-?+-Pyiv^e'('^^-''^*=f/^)'^) dA;, (10) 

? J-m/2 

which, as a function of a and r, only depends upon the parity of rq + ap. To obtain the next 
formula we substitute the estimate r^e^^'' ~ e~^^^ e^'^^^ "^^ (where we take k mod m 

to lie between — m/2 and m/2). Doing this substitution introduces an error, but we postpone 
the error analysis until later. Note that we can substitute this approximation for r^.e"^'= cither 
in (10) or just prior to the integral approximation (9), since the only property of r^e*''^'" that 
we used in the intervening steps is that it is periodic in k with period m. 

9 J-m/2 

We extend the range of integration to all of M, introducing a negligible error. 



{qn) 



J —00 
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To measure the closeness of (6/(27r) top/g, we define a so that (j)/{2'K) = p/q{l + aW), where 
we define W = ^/qe/ (np). Then (0 — 2'Kp/q)qk — 2'KpaWk = 2ay/qek. We change variables 
to a; = ^/qek to obtain 



/oo 
Lii_^ (^(_l)rg+ap^<jg2iax-x2-) ^^3^ 
-00 



The reader may wonder about the apparent asymmetry in equation (13) (when £ = 0), 
e.g. why does q appear but not p, while Z is symmetrical with respect to width and height? 
But when h -\- c ^ a 1 and {nh) / {mc) ^ p/q we have 



W 



qe riT nr ^/2 



np V P'^c V {pqmnbc) ' 

1 _ 1 _ {mnhcfl^ 

log/1^ = nglog r ~ (o + c — a) — ~ (o + c — a)——^ 

a — b c 



(14) 



{qnY (pqmnc/bY^'^. 

The asymmetry between b and c when i > should not be unexpected since we differentiated 
with respect to c rather than b. 

Referring back to the spirals in Figure 4, q counts the number of spirals, the parameter 
A measures the radius of the spirals, ^/e is a measure of how far apart the points are on the 
spiral, and a is a measure of the "spirality" . a is the right parameter against which to plot 
the shape of the spikes, making W a measure if their width. When a — the spirals are 
spokes, when a gets too large (for a given ^/e) the spirals break up into a cloud, by which 
time the integral approximation (9) breaks down. 

4.2 Error analysis 

In the interest of simplicity, we only consider the case when p, q, A, and a are held fixed 
while m ^ oo and n oo. Most of the interesting behavior already shows up in this case. 
It is also quite interesting to ask how much these parameters can vary (e.g. can p and q be 
as large as n^/^?), but we do not pursue that in this article. 

Lemma 1 When we fix p, q, A, and a while m — > oo and n — > oo, we have 
■J^e d d 



{qny 



a a poo 

C---- C- log(Z,0 = -J Lil-^ ((_l)r9+-P^.e2.a.^x-^ ^ ^^^^^ 



provided that the curve (^—iyi+^PA'ie'^^°'^~^^ does not contain the point 1. If the curve does 
contain 1, then the convergence for i — is still valid provided that no multiplicand of Z^^ 
is closer than e~°^^^ to 0, in which case the right-hand side is merely an upper bound. 



Critical Resonance 



Kenyan Sz Wilson 



13 



These integrals are explicitly evaluated in § 5.1. 
Proof: Much of the proof has already been given in § 4.1, what we have left to do is 
justify the approximations that we made in (9) and (11). For this error analysis we do the 
r^e^'^^ ~ e"*^*^ substitution before the integral approximation. With k = j + qu and 
jeZ, + a/2, 

— Lii_^ ^(^_]^jTg27riip/g^g2iay^A;-efc2+0(fca/m2)'j 

using — 2'Kp/q — 27r{p/q)aW — 2a^Je|q. Next we use the fact that for integer ^ > 0, 
j^Lii_i{e^) = Li_£(e^) is a rational function of with a pole ai — 1 but which is 
bounded outside a neighborhood of this pole, 

which is vahd as long as (— l)^e^'"-^^/^Ae^*"V^'=~^'=^ lies outside a neighborhood of 1, and 
the 0{k^ /rri^) error term is much smaller than the radius of this neighborhood. 

It is not hard to see that the q curves (— l^'^e^'^'-'^/^74e^'"'V^'^~^'^^ are bounded away from 
1 if and only if the curve (^—lyi+'^p ^^lax-x avoids 1. 

Adding up the errors 0{k^/m^) over the range \k\ < Q{^/n) gives 0(1). When \k\ ^ 
O(y^), both and its approximation e~'^^ are exponentially decreasing in |A;|, so using 
\jii_i{z) ~ z for small and = e^"^ (1 + 0{k^ /m?')) when k < n^/'\ we sec that doing 
the substitution for n^/^ < k < r?!'^ also introduces 0(1) error, and that the substitution for 
k > r?!'^ gives o(l) error. Upon multiplying by ^/^fe/^qnY, all these errors become o(l). 

For the integral approximation (9), the integrands in (9) are continuous (except at the 

branch cut) as long as the curves (— l)^e^'^*-'^/''Ae^*"v^^''^''^^ avoid the point 1. Moreover 
they converge exponentially fast to when \k\ oo. Therefore they are Riemann summable 
and the error in converting the sums to integrals tends to zero. The error introduced by 
extending the range of integration to the reals is exponentially small. 

What happens when a curve passes through the singularity? When £ = 0, the integral in 
expression (9) for \ogZ„j. converges and is an upper bound for log Z^^r'- the Riemann sum for 
log Z^T- converges to its integral on the complement of a small neighborhood of the singularity, 
and the Riemann sum near the logarithmic singularity has a negligible contribution except 
possibly for the point which is closest to the singularity. If the distance of the closest point 
to the singularity is no smaller than e~°^^\ then the contribution of this point is o{^Jn) and 
so can be ignored. □ 



4.3 The distribution of the number of edges of type 'c' 



By Lemma 1, to first order logZ^j- only depends on the parity of rq + crp, so we define 

1 
2 



^- = \ X] ^o-r and = \ X] ^arZar 



Cr.T (T,T 

Tq + ap odd rq + ap even 



where e^r = — 1 if cr = and r = 0, and e^r = 1 otherwise. We have Z — Z^ + Z^, and from 
Proposition 1, both Z_ and Z+ have only nonnegative coefficients, so they can be interpreted 
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as distributions. Lemma 1 shows that typically one of Z_ or Z+ is exponentially larger than 
the other one, so the distribution of the number of c edges is governed by whichever of Z_ 
or Z+ is dominant. 

We use the method of moments to determine the distribution of the number iVg of type-'c' 
edges. We saw in § 2 how to express the £th moment of Nc about its mean (call it Ci) in 

d d d d 

terms of — c— ■ ■ ■ c— log Z, but in Lemma 1 we evaluated Kg ^-t = c— ■ ■ ■ c— log Z^jt. 

oc oc ' oc oc 

Define Ce^ar to be the same expression, except with the Ki^^ttS replacing the K^s (see also 

(15) below), and similarly define the Ki^±s and the C^,±'s. 
Lemma 2 Under the assumptions of Lemma 1, if the curve 

contain the point 1 and Li_i ((_i) Wp^gg^^^^^-x^-j _^ then C^^^rlC^-Ilr ^{^- 1)"- 
(Eere as usual £\\ — i{£ — 2) • • • (3)(1) when i is odd and i\\ — when I is even.) 

Proof: Prom Lemma 1 we have Kj^^^ — 0{{qny / -^/^e)., and since the above integral is 

nonzero, C2,aT = -?^2.o-t = QHqn)^/ Thus each monomial in the polynomial (1) for 
C^^CTT has magnitude 0{{qnY / {^J q'^e)'^^'^^''^'^) . Recall that q'^e ^ 1. The monomial degree is 
uniquely maximized by the K2^^^ term [i even) or the K^^^rKi^a?'^ term (£ odd). Thus 
when £ is o dd, Ct,^r/C^£r = 0((g3e)i/4) ^ 0, and when f is even, Ct^arjC^I-r tends to the 
coefficient of the monomial K^"^^ in Q^o-t) which by (2) is (£ — 1)(£ — 3) • • • (3)(1). □ 

Thus the C^^o-r's converge to the moments of a Gaussian, but recall that we cannot view 
Zo-T as a distribution since it may have some negative coefficients. Next we show that the 
moments Ci^± of Z± (which are genuine distributions) are close to the corresponding C^^o-t-'s. 
Lemma 3 Suppose £ e = (-l)^i9+^i?' = (_i)^2g+<x2P^ (ai,ri) (a2,r2), and either 

q — —1, or else = +1 hut A>1. Under the assumptions of Lemmas 1 and 2, if\Ki^^^j-^ — 
Ki,,,r,\ « n'/^ then Q,,/C^/,%^ = (^ _ 1)11 + o(l). 
Proof: Viewing Z^t- as a polynomial in c, Z^^t- — 'Y1ii1i,tt c^ we have 

Ci,cTZaT = li,<TT (^ - Ki^arYd. (15) 

As Z^ is the average of So-iTi-^o-in ^-nd Scr^j-^Z^^j.^, we have 



2 

cr,T 



1 ^ ft\ 

(T,T 1 j=0 ^-^^ 

= 2 ^'^'^ \j) '-'3,'^rZcT{Kl,aT - Ki^<;Y~^- 

Since = — 1 or else <; = +1 but A > 1, ^o-in^o-iri and £^2x2^0-2x2 have the same sign, 
so Ki^^ is a convex combination of fTi^o-iri and Ki^cr2T2i and it too can differ by at most 
< ^ 0{Cl^^X) from them. Substituting Cj^^r = ((j - 1)!! + o(l))C(/^^^ we get 
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so Ci^^ is a convex combination of ((£ — 1)!! + o(l))C2^^^^j and ((^ — 1)!! + o(l))C2^^2r2- 



Since C2,o-iti — ^o-in and C2 



K^^a^T^i Lemma 1 and the hypothesis on the A'2's 



£/2 



cm • 



□ 



imphes C2,.,., = C2,.,r, + 0(n5/4). Thus C,,, = ((£ - 1)!! + o(l))C2' 

We have not computed i^i^o-T to the precision that Lemma 3 would appear to suggest 
that we need, but aU we reaUy need is that the two relevant i^i^o-r's are quite close. 

Lemma 4 Suppose ^ = (^—lyfi+'^iP — (^_]^^t2?+o-2p Under the assumptions of Lemma 1, if 



2iax—x 



then \Ki 



Lii(-^^^e' 



I < exp(-0(v/^)). 



dx < — 



Lii(<j^% 



2iax—x 



dx, 



Proof: For expository convenience say that the two cr's for which {^—lyi+'^P = q are 01 
and 10. From Lemma 1, Zqi and Zio dominate —Zqq and Zu by a factor of exp{Q{y/n)), 
and then from Proposition 1, \Zoi/Ziq — 1| < exp(— 9(y^)). Writing Z^t as a polynomial 
in c, Z^T — Yliili,crT d' 1 again from Proposition 1 we have 

|7j,oi - 74,10 1 < -7i,oo + 7j,ii 
|7»,oi - 7j,io ic'\ < «c*(-7,,oo + 7i,ii) < "^^^(-74,00 + 7j,ii d) 



< mn 



|^oi-f^i,oi - ^10-^^1,10 1 < mn(-Zoo + ^n). 

As Zqi and Zio arc exponentially close to each other and exponentially dominate —Zqq and 
Zii, it must be that ^^1,01 and Ki^iq are exponentially close. □ 



Theorem 2 When we fix p, q, A, and a while m — > 00 and n — > 00, we have 

{mnbcY^* 



\ogZ^ 



7rA/2(pg)3/4 



max — 

± 



Lii ( ± dx + 0(1) 



(16) 



and with the exceptions noted below, the number Nc of edges of type 'c ' converges in distri- 
bution to a Gaussian, with mean 



[ranc 



1 3/4 



and variance 



TTV2{pqby/^ 
{mncY^^{pqY^* 



/oo 
Lio(±AV'"=^-^')dx + o(l) 
-00 



/oo 
Li_i (±AV^"^-=^')dx + o(l) 
-00 



(17) 



(18) 



where the choice of± in (17) and (18) is the value that maximizes (16). The exceptions are 

1. When both + and — maximize (16), the distribution of type-'c' edges is a mixture of 
the two Gaussians defined above, provided that exceptions 2 and 3 do not also occur. 
(In the interest of space we omit the proof about the mixture of Gaussians, but we can 
supply it to the interested reader upon request.) 
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2. If for the dominant choice of i: in (16), the curve ±^9e^*"^~^ passes through 1, we 
do not say anything about the distribution of type- 'c ' edges, but the formula for log Z 
remains valid. (We have reason to believe, but have not proved, that this scenario never 
occurs.) 

3. In the event that the integral in (18) evaluates to 0, the formulas are still valid, hut 
we no longer claim that the distribution is a Gaussian. (We believe that this scenario 
never occurs. Lemma 6 in ^ 5.2 rules it out when A<1.) 

(Sec Appendix B for a review and discussion of the parameters A and a.) 
Proof: Immediate from Lemma 1, Lemma 2, Lemma 3 (with the fact that Lemma 5 
shows the dominant choice of ± to be — when A < 1), Lemma 4, the fact that the moments 
are those of a Gaussian random variable, the method of moments, approximation (14) for 
1 / y/ q^e, and qn ~ y/pqmnc/b. 

In the event that one of the curves ±^9e^'"^~^ passes through the point 1, we still obtain 
(16) from Lemma 1 because each such curve has two corresponding Z^t^s, at most one of 
which can have a multiplicand closer than e^"*^^-* to 0. □ 

Theorem 1 follows by plugging into Theorem 2 A = 1 and a = 0, using the explicit 
evaluation of the integrals in § 5.1, and using — Li3/2(— 1) > — Li3/2(1), — Li_i/2(— 1) > 0, 
and Li,(-1) = {2^-" - l)C{i^). 



5 Understanding the resonant spikes 

The resonant spikes, an example of which is shown in Figure 5, exhibit nontrivial behavior. 
We saw already that this behavior is determined by the integrals 

/ Li, (/3e2^°^-^ ) dx, 

J —oo 

with u — 1 governing logZ, u — governing (A^c), and u — —1 governing a'^{N^. We start 
by explicitly evaluating these integrals in § 5.1, and then investigate some of their properties 
in § 5.2, § 5.3, and § 5.4. These subsections may be read in any order. 

/oo 
Li^ [(5e^'"^) dx 
-00 

We assume that /? G C, « G M, and G C, though later we restrict v to integers < 1. For 
convenience we assume a > since the integral is an even function of a. When \ j3\ > 1 there 
is a branch cut (Li, (2;) has a branch cut {2; G R : z > 1}) that may be encountered when 
we vary x, so we need to specify which branch of the poly logarithm we are integrating over. 
Since the principal branch is the one for which Li, (/^e^*"^"^") ^ ±00, we specify 

that the integrand is the principal branch of Li,, even though this may make the integrand 
only piecewise analytic as a function of x as it ranges from —00 to 00. In the event that 
q; = and (5 is real and > 1, so as to ensure continuity in a, we specify that ^e^'"^~^ lies 
above the branch cut for positive x and below the cut for negative x. 
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log Z{j)qf'l^ 



when A"^ 



when A*^ 



a 



when A*^ 




-3 -2 -1 



log Zjpqf I ^ 

— , ,^ ,, when A^ = exp(l) 

1okZ_ loffZ, loffZ_ \ogZ+ logZ- 



logZ_ logZ+ logZ_ logZ, logZ_ 




-2.7397r -2.2377r -1.7357r -1.2337r -.7407r -.29l7r .29l7i .7407r 1.2337r 1.7357r 2.2377i 2.7397r 



t t 



Yj^ when ^4^ = exp 



11 



^ ^ ^ ^ ^ 



-2.7397r -2.2377r -1.7357r -1.2337r -.7407r -.29l7r .29l7r .7407r 1.2337r 1.7357r 2.2377r 2.7397r 



Figure 5: Anatomy of the resonant spikes. In the upper panel we show the curves for 
logZ, (Nc), and a'^{Nc) as a function of a when A = 1. (The parameters A and « are 
reviewed in Appendix B.) When A < 1, Z_ is dominant, and all three curves are analytic 
and unimodal. The situation is quite different when A > 1. In the next panel we show 
the curves for logZ+ and logZ_ as a function of a when A"^ — exp(l) (logZ is the max 
of these two curves). When A > 1 there are singularities in logZ_|_ and \ogZ_ that occur 
when their corresponding spirals cross the singularity, that is, when a = (tt/ v^log Ai)Z for 
logZ+ and when a = {n/y/\og A'?)(Z + 1/2) for logZ_. In the lower two panels we show the 
curves for (Nc) and a^{Nc). These curves were computed using (17) and (18) as explicitly 
evaluated in (21) with whichever of Z+ or Z_ is significant. The "crossover points", where 
log Z_ and log Z+ alternate in significance, are in a sense a phase transition within a phase 
transition. At each crossover point, the curve for logZ is continuous but nonanalytic, the 
curve for (Nc) is discontinuous, and the curve for a'^{Nc) has a delta function. When a spiral 
hits the singularity, the corresponding Z± appears to be the insignificant one. The large-a 
asymptotics of the curves for logZ, (Ac), and cr^(A^c) are given in Theorem 3. 
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We z = x-\-iy (x, e R) and integrate instead within the complex plane. Rather than 
integrate \Ay (/3e^*"^~'^ ) along the real axis 'isz = 0, it is more convenient to the integrate 
along the line — ol. When we deform the contour of integration and set z = x + ia we 
get 

/oo foo terms from the singularities and 

Li^ (/3e2*°^-^') dx= \ Li^ (/^e-^'e"^') dx + branch cuts of Li^ (/^e^^"^"^') in 
"°° °° the complex 2;-plane. 

In particular if < 1 there are no singularities or branch cuts encountered when de- 
forming the contour of integration, so there are no additional terms. We shall evaluate the 
main term in first. For now assume that |/?|e~° < 1 so that the series expansion of the 
poly logarithm is absolutely convergent. This enables us to write 

°° ( o^—a^Xn roo 



n=l * 

= V7rLi^+i/2 (/3e~°'). 
The singularities of hi^, (/3e^*"^~^^) in the complex ^-plane occur when 

P^2iaz-z'' = 1 

2iQ;2; — z^ — — log /3 — 27ri/c 



z — ia± %\/ — log/3 — I'Kik 
where k is an integer. Since we are moving the contour between = and Qz = a (see 



Figure 6), the relevant singularities arc of the form ia — i\/ — log/5 — I'nik where the 
principal square root is taken. (If arg/3 = 0, < log|/3|, and /c = then both roots are 
relevant.) The branch cuts of the integrand occur where /3e^'"^~^ is real and > 1. With 
z — X ^iy and y e R these are 

2ax — 2xy = — arg f] — 27rk —2ay — x'^ + y'^ > — log \P\ 

127rA; + arg/5 2 2 ^ 2 , 1^1 

x = - [y-a) -X >a -log|/3| 

2 y-a 

To identify if and where the kth cut intersects the line '^sz = we set y = and find that 
there is an intersection at 

2nk + arg /3 
2a 

provided 

(27r/c + arg/?)2 < Aa"^ log\(3\. 
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contour for 

Li^ (^e-"'e-^') 





Figure 6: Singularities, branch cuts, and contours of integration for Ia^ (/3e^*"^~^^) . When 
pushing the contour from = a to $52; = 0, the contour must deform to travel around the 
singularities and branch cuts, so these contribute to the integral. 



For the moment suppose that |/9|e~" < 1 so that there are no singularities or branch cuts 
of the integrand on the line = a. Push this contour downwards toward the line = 0, 
except that when a branch cut is encountered, the contour must go around the branch cut, 
and travels along the cut upwards on its left side in downwards on its right side (Figure 6). 
The integrand increases by (27ri/r(i/))(log(/3e^'"^~^ )y~^ when going from the left side of 
the cut to the right side, so the net contribution of this cut to the integral is 



27ri 



(log(/3e 



2iaz- 



')Y-Uz 



c 



(plus another term from the singularity) where the contour C is a branch cut travelling 
from the branch point to the point where the branch cut intersects the line '^z — 0. For 
the kih. branch cut the endpoints of integration are ia — i\l — log /9 — 27ri/c and —{27rk + 
aigP)/{2a). When u — 1 the singularity does not contribute and this integral is easy to 
evaluate and we get 



27ri 



+ i\/ cP- — log /3 — 2ix%k 



27rk + arg (3 
2a 



Since this term gets added to the integral over ^z — 0, when we evaluate the integral over 
^z — we subtract this term from the integral over ^z — a. Thus 



Lii {(3. 



27rz J2 



-la 



+ i\/ o? — log/3 — 2'K%k 



{2iTk+axgl3f<4a^ log | 



27r/,- + arg (3 
2a 



(19) 
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There is a term in the summation for each time the spiral /5e^'"^~^ encloses the singularity 
at 1. This formula is valid when |/5|e^" < 1, and we would like to show that it is valid 
(for real ji and a) without this restriction. We can do this via analytic continuation in 
(complex) /? and a. While both sides of equation (19) arc analytic in a (except where the 
spiral ^e^*"^~^ crosses the singularity), neither side of (19) is analytic in /3. The arg/3 terms 
on the right-hand side are nonanalytic, and the left-hand side is nonanalytic due to the 
branch cut in Lii(^) — as arg (5 is varied, portions of the integrand cross the branch cut. To 
remedy this problem, we deform the branch cut as arg (5 changes, so that the cut consists of 
the arc {2; G C : |^| = 1 and < argz < arg/3 or > arg 2; > arg/3} together with the ray 
{2; e C : > 1 and arg 2; = arg/3}. So long as the spiral avoids the singularity, for any 
value of X, the integrand does not cross the moving branch as /3 and a are changed. This 
re-interpreted integral (with the moving branch cut) is then analytic except where the spiral 
l^^iax-x CTosses the singularity, and the arg /? terms on the right-hand side become constant. 
The \/7rLi3/2(/5e~" ) term on the right-hand-side of (19) is nonanalytic when /3e~" = 1, but 
this nonanalyticity is cancelled by the nonanalyticity in the A; = term in the summation 
(by (24) in Appendix A, which reviews polylogarithms) . Thus the right-hand-side of the 
re-interpreted (19) is also analytic whenever the spiral ^e^*"^~^ avoids the singularity. For 
a given /? and a, we may continuously decrease \i5\ and increase a until |/9|e~° < 1 (where 
we already know that (19) is valid), while keeping the spiral from crossing the singularity. 
Thus (19) is valid for all complex /3 and a for which the spiral avoids the singularity. As 
both sides of (19) are continuous, (19) is also valid when the spiral contains the singularity. 
We revert to the principal branch cut interpretation of the integral, restoring the arg /3 terms 
on the right-hand side of (19). 

When /3 is real, upon summing over A;, the terms (27rA;-|- arg/3)/(2Q;) cancel: when /3 < 
the /cth term cancels the (—1 — A;)th term, and when /3 > the A;th term cancels the — /cth 
term, while the 0th term does not contribute. For real /3 we get 

/ Lii (/3e2*"^-^') (ix= V7^Li3/2 (/3e-"') -27r ^ [a - ^ -\og^ - 2mk\ . (20) 

(27rfe+arg/3)2<4a2 log |/3| 



When /3e " > 1 the k — ^ "branch-cut term" exactly cancels the imaginary part of the 

(A 

9/3 



Li3/2 term. Applying multiple times to (20) we find for real /3 and nonpositive integer v 



/oo 
Li,(/3e2^"^-"')(ix = V^Li,+i/2(/3e-"')- V^r(i-z/) ^[a^ _ log/3 + 27riA;]-^/2+._ ^31) 
-00 u^nj 



(27rA;+arg^)2<4a2 log | 



5.2 Inequalities 

To use Theorem 2 we need to know which of Z± is significant. We already saw that they 
can alternate in significance (as a is varied) when A > \. We claim that when A < 1, or 
else A>\ but a = 0, that Z_ is the significant one. 

Lemma 5 //O < /3 < 1, or else /? > 1 hut a = 0, then 

/oo /-oo 
Lil (/3e2^"^-^') dx<- Lii ( - /^e^^"^-^') dx. (22) 
-00 J —00 
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Proof: Under either of the hypotheses, equation (20) simphfies to 

/CXD /"OO 
Lii ( ± /3e2'"^-"') dx = ( ± /?e-"') = / Lii ( ± /5e-"'-"') dx, 

-oo J —oo 

but 

-3?Lii(/3e-"'-^')= 3?log(l - /^e""'-^') < 3?log(l + /3e-"'-=^') = -3?Lii(- /3e-"'-^'). □ 

To show that the distribution of N^. is a Gaussian, we need to know that the integral 
expression in (18) (for the dominant choice of ±) is nonzero. The integral clearly cannot 
be negative, because it has an interpretation in terms of variance, but a priori it could be 
zero, in which case the o(l) error term would control the variance, and we would be unable 
to characterize the distribution of Nc- Ideally we would like to show that it always positive, 
but we only do this for >1 < 1, or else A > 1 but a — 0. Recall that under these conditions 
Z_ is dominant. 

Lemma 6 7/0 < /3 < 1, or else /3 > 1 but a — 0, then 

/oo 
Li_i ( - /3e'^"^-^') dx > 0. 
-oo 

Proof: As above, under either of the hypotheses, equation (21) simphfies to 

/oo /"OO 
Li_i ( - /?e2'"=^-^') dx = A/^Li_i/2 ( - /9e-"') = / Li_i ( - /^e""'-^') dx, 
-oo J —oo 

but - Li_i ( - /3e-°'-^') = /3e-°^'-^V(l + /Se""'-^')^ > 0. □ 

Lemma 7 The curve for logZ attains its (unique) maximum value when a — 0. 

Proof: When < /3 we have log(l + fJe"^^) > 3?log(l ± /3e^*"^~^^), with strict inequality 
when ±6^^"=^ 1, so - Lii ( ^ /^e^^"^-^') dx < - J^^ Lii ( - /^e"^") dx, with strict 
inequality unless both a — and ^ — —. □ 

In contrast, the curve for (Nc) does not always attain its maximum value when a — 0. 
5.3 Asymptotics for large a 

In Figure 5 the curves for logZ and (Nc) appear to asymptote out to a positive constant, 
while the curve for a^(A^c) decays to 0. The following theorem gives the large-o; asymptotics 
of these curves which are given by the integrals in (16), (17), and (18) in Theorem 2. 

Theorem 3 For positive (3 and real a, 



3(l0g/3)=^/2 /3>1 



/OO I 4 

Lii (±/3e^*"^-^') dx = 
-oo ^ ^ (0 /9<1 

lim-rL.o(±/^e— 

H-oo ^ ^ jo (3<1 

(3<l 



hm -a' / Li_i ( ± /Je'^"^-^ ) dx 

\a\-^oo 
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Proof: These integrals are evaluated in (20) and (21). When \a\ — > oo, the main terms in 
(20) and (21), — y7rLii,+i/2(±/9e~" ) for t/ = 1,0, —1, become negligible. Thus we only have 
to compute the "branch-cut terms" when (5 > \. For (20) these terms can be approximated 
for large a as follows. Let 7 = log/5 — o? . The identity 

■\/x-\-iy + yjx-iy = \ll{x + sj + \p) 
is easily verified by squaring both sides, and allows us to rewrite the sum in (20) as 

^ 27r {\a\ - ^(-7 + A/? + 47r2A;2)/2') 

where k is integer for positive (for Z_|_) or half-intcgcr for negative (5 (for Z_\ 

Suppose |q;| S> v^log Then /c, which ranges up to (|Q;|/7r)-\/log/9, is much smaller than 
I7I = (1 + o(1))q;2, and we have 



7 + v/7T4^)/2 = (1 + Vl + 47r2A;V72)/2 

= + + o(l))F/72 



v/«2_iog^(l + (^2/2 + o(l))fcV7') 

|a|(l - (1/2 + 0(1)) log/3/a2)(i + {j^ + o(l))A;Va^) 

|a|-i±4^(log/?-7r2^Va2). 



In particular 



27r (^|«| - ^(-7+A/72 + 47r2A;2)/2^ = !L±^(log/3 - Tr'/tVa") 

and we sum this over /c's for which the quantity is positive. If \a\ ^ 1 / \/\o'g /3 then we can 
approximate this sum with an integral, and the integral is 

7r + o(l)2 2\a 



log/3 x^yb^=(l + o(l))-(log/3)^/^ 

|a| O TT o 

These are the asymptotics for v = \. For integer i/ < the range of h is the same as it 
is for i/ = 1, i.e. up to about \a\\J\o^(i jix . Thus we need to evaluate 

0Fr(l/2 -v) J] [a^ - log/3 + 27riA;]-^/'+^ 

|fc|<|a|VI^/7r 

For q; ^ f3, A; <C so that 

[a^ - log/3 + 27riA;] + [a' - log /3 - 27riA;]-^/2+'' ?a 2|q;|2''-^ 
and if a ^ 1/ ^/\ogP the summation is approximately an integral which is asymptotically 



^r(l/2 - u)a"'^/h^. 
Taking u — and 1/ — —1 give the desired results for the (Nc) and a'^{Nc) integrals. □ 
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5.4 Crossover locations 

\i A <\ then Z_ always exceeds Z+. But for larger A there are crossover values for a at 
which Z_ and Zjy. alternate in significance. Since Z_ and count configurations in different 
Z2-homology classes (see (4)), we conchide that the crossover values arc the places where 
the typical homology class of a lattice path changes. So there is a phase transition at these 
points: the topology of a typical configuration changes. 

The crossover values a for a fixed ^ — satisfy an imphcit (and transcendental) equa- 
tion. Instead of solving these equations directly we can derive an analytic expression for the 
crossover a's as a function of 7 = log/? — ct^, and then given 7 and the crossover a we can 
calculate the corresponding /9. In this way we can parametrically plot these critical pairs 
(/3, a) as a function of 7. For example, for the 0th crossover we have 

^l/2^3/2^_ + 0(1) = -v/^Li3/2(-/?e-«') 

^1/2^3/2^^ + 0(1) = -V^Li3/2(/9e-"') + 27r(a - V«'-log/?) 



from which we can solve 



Li3/2(e'^) - Li3/2(-e'^) 

" = 2A + 

Similarly, for the rth crossover we have 

I Li3/2(e^) - Li3/2(-e^) ^ ^ , ——, 

Next let us approximate these crossovers for r fixed and 7 large. For large 7 we can 
substitute z/ = 3/2 into the asymptotic series expansions (25) and (26) for Li,^(±e^) to write 

- Li3/2(-e-) « + \^'l'h'l' + ^/''h"' + ■■■ 



so that 



Li3/2(e^) ^ —y^V^' - 1^1^ ^\7^'l^ h"^ + ^/'^h"'" + 



r 

(-1)''q; ^7~^/^ + J^(-l)^ V'-7 + ^TTi + y-7 - /cTTi] 



A;=l 
r 



fc=i 



77-'/' + $](-l)V2(-7 + 7(l + ^%V(V))) 



4 

fc=i 

r 

^7-^^' + E(-i)'^/^^ 

fc=l 

f 7-^/^ + $:(-i)^wy/^ 

fc=i 

(-l)^(r/2 + l/4)7r7-i/^ 
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Since a. is small, 7 log/3, and so the rth crossover occurs at a (r/2 + l/4)7r/v/log/3. 

By comparison the nonanalyticities in the curve for Z_ occur exactly at a = vr/ v^log /3(Z+ 
1/2), and the nonanalyticities in the curve for Z+ occur exactly at o; = ixj^/Xog 

6 Open problems 

Our analysis for of the lattice paths at the critical point is geared to regions whose aspect 
ratio is close to a simple rational number, but we do not know how to treat "irrational 
domains" which do not have a simple rational approximation. Consider for instance the case 
a = l, 6 = c= l/2 when the side lengths m and n are successive Fibonacci numbers. Then 
the aspect ratio of the region is very close to the Golden ratio, which is not well approximated 
by simple rationals. In this case we believe that the partition function Z is 0(1) (about 2.1), 
so that with 0(1) probabihty there are no lattice paths. In the event that there are lattice 
paths, we believe that they connect up into 0(1) loops each of length 0(n^/^). We have also 
found empirically that as one varies the aspect ratio of large regions, the smallest value that 
the partition function Z takes on is very close to 2. We do not know how to prove any of 
these conjectures. 

In Figure 3 we plotted logZ and (A^c) versus njm when a — \^h — c— 1/2. As predicted 
there are spikes in logZ and (iVc) when the aspect ratio njm is near a simple rational 
p/g. But there also appear to be flanking secondary spikes in (A^c) when n/m ~ p/g but 
|n/m — p/g| ^ Xj \fn. We do not understand this phenomenon. 

When a = 1, b = c = l/2 + 0(l/n) and the aspect ratio is nearly a simple rational 
p/q, the partition function and edge density as a function of a are nonanalytic at certain 
points. We know that the homology type of the strands changes at these nonanalyticities. 
We conjecture that when a = each loop winds around exactly p times horizontally and q 
times vertically. For large enough values of a, when one follows a loop for mp + nq steps, 
one does not return to the starting point. We conjecture that the number of nonanalyticities 
between a and determines how many strands away one ends up after following a loop for 
mp + nq steps. 

Despite our explicit formulas for \ogZ± and for the expected value and variance of the 
number of edges, there are some basic properties about these functions that we have not been 
able to derive. For example, we conjecture that when our formula for log Z_ is nonanalytic, 
our formula for log gives a strictly larger value, and vice versa, and that the formula for 
a'^{Nc) is strictly positive. We also conjecture that ior A > 1 our formula for logZ as a 
function of a always exceeds its limiting value as o; 00. These conjectures seem fairly 
evident from the graphs in Figure 5, but it is not obvious how to prove them. One would also 
like to determine for fixed A > 1 the maximum and minimum values given by our formula 
for {Nc), since both of these can be different than the limiting a — > 00 value. 
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A Polylogarithms 

The poly logarithm function Li^{z) is defined by 

n=l 
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for l^l < 1 and by analytic continuation elsewhere. We may for instance write Lii(z) = 
— log(l — z), Lio(-2;) = -2/(1 — and Li_i(2;) — zjiX — zf. The polylogarithm function has 
the convenient property that 



Appell's integral expression 



L'^W = rM/. — '''' 

is valid for 3?(i/) > and z ^ [1, oo), and defines the principal branch of the polylogarithm. 

The polylogarithm has an interesting Riemann surface. When one crosses the branch cut 
[1, oo) in the positive direction, the polylogarithm increases by ^(log,2)'^-^ (Both this and 
the defining scries expansion for Li,^ arc readily derived from Appell's integral expression.) 
For nonpositive integer u this quantity is 0, consistent with the fact that Lii, is a rational 
function for these values of u. Unless u is an integer < 1, the Y^{^ogzy~^ term creates a 
second branch point at 2; = off the principal branch. The function Li^{z) is analytic in 
both z and u, except for a singularity at z = 1 (and z = oS the principal branch). For 
further background see Bateman and Erdelyi et al. [1, Chapt. 1 §11], Truesdell [28], Dingle 
[10, 9], and Lewin [16]. 

When z is on the principal branch and near 1, the series expansion 

U^{z) = r(l - log zf-' + E - n)^^^ (24) 

n=0 

was given by Lindelof [17, pp. 138-141] (derivations are also given in [28] and [5]), and is 
absolutely convergent when |log^| < 27r. 

We also use the asymptotic series expansions for lAy{z) when \z\ 00. For large positive 
X these are 



°° (1/2^''"^ - 1) C(2A;) 
Li^(-x) = - cos(7ri/) Li^(-l/a;) + 2 J]] V{v + l-2k) — ^^^^^ 



-2k 



k=0 

and 



(25) 



Li.(a;) = - cos(7r^) Li.(l/x) ± ^^^^ + 2 g r(z. + 1 - 2fc) ^^"g"^^'^"''" ^^^^ 

Aside from the — cos(7rz/) Li,y(±l/a;) terms, the asymptotic series expansions (25) and (26) 
were derived by Sommerfeld and Clunie respectively. Rhodes first derived (25) for the case of 
integer u, and for these u the expansions (25) and (26) have only finitely many nonzero terms. 
For noninteger u the expansions (25) and (26) diverge, and the — cos(7ri/) Li,^(±l/a;) term is 
dominated by each term in the divergent series expansion, so one may wonder what role it 
plays. Dingle [10, eqn 17] [9, eqn 11] showed how to make practical computational use of 
these series by truncating the series after finitely many terms and providing convergent series 
expansions for the remainder when the — cos(7rz/) Li,^(±l/a;) term is present. For positive 
X > 1 we are evaluating Lii,(x) on the branch cut of the principal branch; the ± sign is 
positive when the branch cut is just above the real axis, which is the usual convention. See 
also Pickard [21, eqn 3.5] for the asymptotics for large complex values of z. 
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B Notation 

a, 6, c: weights of edges in the three different directions, a is weight for vertex not being in 
a loop, h is weight for horizontal step, and c is weight for vertical step. 
m: horizontal length of torus 
n: vertical length of torus 

Simphfying assumptions: h < a, h/ a = 6(1), 1 — 6/a = 0(1), n = 6(m) 
p/q: rational approximation of nb/{mc), gcd(p, g) = 1. Intuitively p is the number of 
horizontal windings of loops, and q is the number of vertical windings of loops, but this 
intuition is not quite accurate if "ratcheting" takes place. 

is a measure of close the rational approximation p/q nb/ (mc) needs to be. For fixed p 
and q, W = Q{l/^/n). More precisely, 

W = ^Jqej {np) = ^2qnab/ {pm{a — b)) ^ a/2/ (pmc) ~ y^2/{qnb) 

a is measure of the error in the approximation p/qK, nb/ {mc) on the scale of W: 

nbq nbq 

mcp m[a — b)p 

A is a measure of how close the weights a, 6, and c are to the phase transition a = b + c. 
When a = 6 + c, and A is sensitive to perturbations on the order of 9(l/n). 

appears in many of our formulas, and is approximately symmetric in the parameters: 



log^'«(6 + c-a),/^^ 



27rn6 
m{a—b) 



= 

= -e^^^ Ok = 2nk/m; k&Z^ + a/2 for Z„ 
r,^ = {a-bY\a + bzk\''' = e-^k^+oi^VrnS) 

4>k = arg(a + 6;2fc)~" = 4>k + 0{k^ /m^) 



